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By combining a squeezed propagating microwave field and an unsqueezed vacuum field on a hy¬ 
brid (microwave beam-splitter), we generate entanglement between the two output modes. We verify 
that we have generated entangled states by making independent and efficient single-quadrature mea¬ 
surements of the two output modes. We observe the entanglement witness Ew = —0.2631q‘q 36 and 
the negativity N = 0.0824+qoo O4 with measurement efficiencies at least 26 ±0.1% and 41 ±0.2% for 
channel 1 and 2 respectively. These measurements show that the output two-mode state violates 
the separability criterion and therefore demonstrate entanglement. This shared entanglement be¬ 
tween propagating microwaves provides an important resource for building quantum networks with 
superconducting microwave systems. 
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When two parties share entanglement many powerful 
quantum communication protocols are available to them. 
For example, they may communicate with security guar¬ 
anteed by physical laws, they may encode data more 
densely than classical bounds, and one party can transfer 
a quantum state to another by transmitting only classical 
information, a protocol known as teleportation [ 1 ]. Fur¬ 
thermore, teleportation can be extended to realize error 
correction schemes [2j. Shared entanglement has been 
a powerful and popular tool for long distance quantum 
communications. A second application for shared entan¬ 
glement occurs in a general quantum information proces¬ 
sor that is structured as a distributed machine compris¬ 
ing many well-isolated copies of a high-fidelity quantum 
register p, 0 • To perform quantum computation, these 
registers must then share entanglement. 

For microwave superconducting qubit circuits, the 
quantum registers that have the longest coherence time 
are built from centimeter-sized microwave cavities con¬ 
taining a few qubits (B). Propagating microwave modes 
are the media being developed to establish and exploit en¬ 
tanglement among such registers. Consequently, entan¬ 
glement between physically distinct itinerant microwave 
modes is an important resource and has been created 
and verified in recent experiments M- Although it is 
possible to verify the presence of entanglement with low 
efficiency measurements 0. 10. o. to perform a proto¬ 
col that exploits shared entanglement between the sender 
and the receiver, such as teleportation or error correction, 
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having higher detection efficiency improves the fidelity 
of the process Bill. Furthermore, the high-efficiency 
measurements of the two parties should have indepen¬ 
dent measurement bases to fully characterize a two-mode 
state. 

For propagating microwave modes the quantities that 
can be measured with the highest efficiency are quadra¬ 
ture amplitudes X and Y, i.e. the cosine and sine com¬ 
ponents of the field relative to some phase reference. The 
two quadrature amplitudes are canonically conjugate ob¬ 
servables; thus, one quadrature can in principle be mea¬ 
sured without added noise, but not both. By adjusting 
the reference phase, one can measure a linear combina¬ 
tion of X and Y. We treat X and Y as random variables 
corresponding to potential outcomes of quadrature mea¬ 
surements. By making repeated measurements on many 
copies of the same two-mode state and adjusting their 
phase references over all possible values, one can fully 
characterize the two-mode state. 

In this paper, we report the generation of entangle¬ 
ment between two spatially separate itinerant microwave 
modes by combining a quadrature squeezed state and a 
vacuum state on a microwave hybrid (Fig. [Tj) . A sim¬ 
ilar method to create entanglement is commonly used 
in optical experiments, such as 0. To verify en¬ 
tanglement, we make efficient, single-quadrature mea¬ 
surements simultaneously on the two separate modes 
with independent control of phases 0% and O 2 . Specif¬ 
ically, we measure quadrature amplitudes of the two 
modes W\(Q\) = X\ cos(#i) ± Yi sin(#i) and = 

X2 cos(#2) + Y2 sin(6 > 2) over all pairs (#1, 62). Because the 
two-mode state in this experiment appears to be Gaus¬ 
sian, we characterize it using its quadrature covariance 
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FIG. 1: (color online) The single squeezer model of the ex¬ 
periment. The squeezer (SQ) prepares a squeezed state with 
squeezing parameter s , where the variance of the squeezed 
quadrature is l/(2s). The squeezed state (red ellipse) and 
the unsqueezed input (green circle) are combined on a quadra¬ 
ture hybrid (HY) to generate entangled modes. The hybrid 
has a power transmission coefficient t and a power coupling 
coefficient 1 — t. The two output modes (orange ellipses) of 
the hybrid propagate onto two physically separate transmis¬ 
sion lines and are fed to the two-channel measurement appa¬ 
ratus to measure quadrature amplitudes Wi(6\) and ^ 2 (^ 2 ). 
The measurement apparatus consists of two single-quadrature 
measurement chains (QM1 and QM2), where each QM em¬ 
ploys a VER as the first amplifier. All sources of loss (in¬ 
cluding loss inside the SQ) and measurement inefficiencies 
are modeled by introducing two fictitious beam splitters with 
power transmission coefficients 771 and 772 , respectively. The 
two squeezed states arrive at the two VERs with fixed but 
uncontrolled phase shifts. We mathematically adjust the ref¬ 
erence phases to align the squeezed states with X\ and X 2 as 
illustrated. 


matrix defined as IY- = (1/2)(ZiZj + ZjZA — (ZAiZj). 
where G {X 1 ,Y 1 ,X 2 ,Y 2 } 

The covariance matrix reveals the correlations between 
the two modes and can be used to demonstrate entangle¬ 
ment. Specifically, one can prove entanglement by ob¬ 
serving violation of the inequality jlbl llq 


R(0i, 62 , a ) = Var 
Var 
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for any nonzero real number a , which accounts for any 
unbalance between the two channels. (We use the con¬ 
vention that for vacuum states Var(FF(0)) = 1/2 for all 
0.) One may optimize the observed violation over phase 
rotations (which cannot change entanglement) and a by 
computing the entanglement witness 

_E\v = min 

0i,0 2 ,a 

so Ew < 0 is evidence of entanglement. For the choice 
of a = 1, we also report Aepr = (l/2)ii(0i, # 2 ,1) which 
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gives evidence of entanglement when Aepr < 1. Al¬ 
though the entanglement witness can detect the pres¬ 
ence of entanglement, it does not measure the amount of 
entanglement. Thus, we quantify the entanglement be¬ 
tween the two measured microwave modes with the neg¬ 
ativity A, which can also be calculated from the covari¬ 
ance matrix [lf| [20|. For bipartite systems, the negativ¬ 
ity is a lower bound on the number of entangled Hilbert 
space dimensions (the Schmidt number) [2lj . 

As is evident from Eq. two vacuum modes are 
not entangled. We must use some device to actively 
transform the vacuum. The Josephson Parametric Am¬ 
plifier (JPA) is the critical piece of technology that al¬ 
lows us to generate squeezed states and to perform effi¬ 
cient quadrature measurements at microwave frequencies 
[22l |. It is a microwave phase-sensitive amplifier that is 
built from an electrically nonlinear microwave resonant 
circuit and that derives its gain from a pump tone excit¬ 
ing the circuit. The JPA amplifies a specific quadrature 
of the input state noiselessly while squeezing the conju¬ 
gate quadrature where the amplified quadrature 

is selected by the relative phase between the pump tone 
and the input state. In this experiment, one JPA de¬ 
noted SQ is used to transform the vacuum to a squeezed 
state. Two more JPAs, denoted VER1 and VER2, act as 
single-quadrature preamplifiers for two microwave mea¬ 
surement chains. By adjusting the pump phases of VER1 
and VER2 separately, we independently control the bases 
of our two single-quadrature measurements. 

In this experiment, we integrate SQ with a hybrid on 
a single chip forming the entangler circuit (25[ (see sup¬ 
plementary for the layout and an image of the entangler) 
to minimize loss in the entanglement generation process. 
Feeding the two inputs of the entangler are two vacuum 
states emitted from two 50 Q terminations thermally an¬ 
chored to a cryostat. SQ is pumped at / s = 6.327 GHz 
and generates a squeezed state with approximate direct 
power gain G s = 3.1 dB and bandwidth B s = 8.5 MHz. 
The squeezed state is displaced in phase space by the 
pump amplitude. This displaced squeezed state and the 
vacuum state interfere in the hybrid, creating entangle¬ 
ment in the two output modes of the hybrid. In order to 
ensure that the following VERs are not saturated, we null 
the SQ pump tone at the input of the VERs with a weakly 
coupled coherent field. The SQ’s gain is chosen as a com¬ 
promise between larger observed squeezing and simpler 
operation of the experiment. SQ gain greater than 5 dB 
will not substantially reduce the measured variance be¬ 
cause of the measurement inefficiency, but the associated 
large pump amplitude will be more difficult to null at the 
VERs’ inputs. 

The two-mode state is then measured by our two- 
channel measurement apparatus. The two output modes 
of the entangler propagate in two separate coaxial cables 
which feed the input of VER1 or VER2, located about 
10 cm apart, forming the first stages of the amplification 
of the two quadrature measurements. The VERs are op¬ 
erated with approximate direct power gain G v = 22 dB 
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and bandwidth E> v = 2.3 MHz, and their outputs are 
further amplified by conventional microwave amplifiers. 
The VERs’ gains are chosen to be large enough to over¬ 
whelm the added noise of the following HEMT amplifiers 
(20 dB) but not larger as that would reduce the measure¬ 
ment bandwidth. Each of the three JPAs’ gains is esti¬ 
mated by measuring its response to a small input tone, 
but with all of other JPAs turned off. Because the JPA 
gains may change by about 2% when all three JPAs are 
turned on, these are only estimates of the gains during 
the entanglement generation and verification. 

Finally, the amplified microwave signals of the two 
channels are mixed down for digitizing with copies of the 
VERs’ pump tones serving as the mixers’ local oscillators. 
The mixers’ intermediate frequency outputs are filtered 
with a 1.9 MHz low pass filter and sampled at 10 MHz, 
yielding new measurements of Wi(0\) and HA(# 2 ) every 
100 ns. The phases 0\ and 62 can be independently ad¬ 
justed relative to each other and relative to the squeezed 
quadrature of SQ. In practice, to adjust 0\ and # 2 , we set 
the pump frequencies of VER1 and VER2 to be 1 kHz 
and 50 kHz above the SQ pump frequency respectively. 
In 1 ms, we acquire 10,000 samples covering the full range 
of both 6 \ and # 2 . We acquire data for 1 s, yielding 1,000 
independent realizations of (Wi, W 2 ) for each pair of (# 1 , 
# 2 )- 

To calibrate our measured variances in units of the vac¬ 
uum, we inject states of known noise into the entangler, 
but we bypass SQ by turning off its pump and operat¬ 
ing it as a noiseless G s = 1 amplifier. States of known 
noise are created by varying the temperature of the cryo¬ 
stat and therefore the 50 terminations that feed the 
entangler. We then measure the variances at the outputs 
of the measurement chains while adjusting the cryostat 
temperature, thus calibrating the measurement chains. 
From this procedure, we also determine that the 50 H 
terminations equilibrate with the cryostat for any tem¬ 
perature above 25 mK (see Appendix El). Because the 
entanglement generation is performed with the cryostat 
below 25 mK, the input variance is indistinguishable from 
vacuum in our experiments. 

In order to interpret the measurements, we construct 
an analytic single squeezer model of the experiment rep¬ 
resented by Fig.[lJ All sources of loss and measured inef¬ 
ficiencies are absorbed into the parameters 771 and 772 . We 
extract parameters by fitting the model to the measured 
variances yielding s = 5.41 ± 0.03, 771 1 = 0.130 =b 0.001, 
and 772(1 — t) = 0.202 ±0.001. Moreover, we are able to 
extract from the model small changes in the VERs’ gains 
when bypassing SQ (see Appendix [A]) . We conclude that 
VERl’s and VER2’s power gains change by linear fac¬ 
tors of approximately —1.70 ± 0.07% and 2.04 ± 0.08% 
respectively when the SQ pump is turned on. 

To test the independence between the two measure¬ 
ment channels, we first examined the quadrature vari¬ 
ances measured separately at the outputs of the two mea¬ 
surement chains. From the measured quadrature vari¬ 
ance of each mode, Var(FFi) and Var(VE 2 ) [Fig.[2fa) and 



FIG. 2: (color online) Separate and joint variances. Shown 
are intensity plots of the measured variances of (a) Wi(0i, # 2 ), 
of (b) W 2 (6>i,6> 2 ), and of (c) \ [Wi(0i, 62 ) ± W 2 (6>i, 0 2 )] cal¬ 
ibrated in units of the vacuum versus the two quadrature 
phases 6 1 and # 2 . (d) An expectation of (c) predicted by the 
single squeezer model represented by Fig. [T| (e) The vari¬ 

ances along the corresponding annotated lines in (c) and (d) 
are plotted versus quadrature phase 0\ of channel 1. The red 
squares are the sum of the green-circle line and the magenta- 
diamond line, where the green-circle line is shifted by ^ in 6 \. 
The arrow indicates the observed value of Aepr < 1. 


(b)], we observe an approximate minimum variance 15% 
below vacuum fluctuation, i.e. 15% squeezing below vac¬ 
uum. Furthermore, Var(Wi), only depends on the mea¬ 
surement phase 6 \ and is independent of the measure¬ 
ment phase O 2 ] likewise Var(W 2 ) depends only on # 2 . 
These plots demonstrate that the two VERs are unaf¬ 
fected by the phases of the other’s pump indicating that 
the two channels are well decoupled. Finally, the 7r rather 
than 27r periodicity of Var(VFi) [Var(W 2 )] as a function 
of 6 \ { 62 ) shows that the SQ pump is successfully nulled 
at the input of the VERs. 

By making joint measurements of the two output 
modes, we detect the correlations between them and 
reveal that they are entangled. In Fig. [2fc) we plot 
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FIG. 3: (color online) Covariance matrix of the two-mode 
state, (a) The covariance matrix calculated from the single 
squeezer model parameter extracted by joint fitting of mea¬ 
sured variances, (b) The 10 independent elements of the co- 
variance matrix are shown for the same data as (a). 


the measured joint variance (l/2)Var(Wi ± W 2 ), and 
in Fig. [2fd) we show the expected joint variance pre¬ 
dicted by the single squeezer model in Fig. [T] Because 
the squeezing we observed from separate measurements 
is diluted with vacuum (Fig. [1]), we anticipate that the 
joint measurements will show more squeezing than the 
separate measurements. Indeed, (l/2)Var(Wi ± W 2 ) 
has an approximate minimum variance 25% below vac¬ 
uum fluctuation. The hybrid generates a two-mode en¬ 
tangled state, distributing the squeezing present in the 
input squeezed state into the two output modes. By 
mathematically inverting the hybrid’s action, the two 
input states can be reconstructed from the joint mea¬ 
surements. For example, only the squeezed input con¬ 
tributes to the variance measured along the magenta- 
diamond line in Fig.[2fc); likewise, only the vacuum input 
contributes to the variance measured along the green- 
circle line. One method of estimating Aepr is to sum 
the variances measured along these two lines [Fig. [SJe)]. 
To see this, note that the second term in the expression 
for Aepr can be written as (l/2)Var [Wi(0[) ± ^ 2 (^ 2 )], 
where 0[ = 6^1 — tt/ 2 and 0' 2 = #2 ± tt/ 2. Thus, if the first 
term in Aepr is evaluated at ( 6 1 ,^ 2 ) on the magenta- 
diamond line, then the second term must be evaluated at 
a corresponding point ( 0 ^ 0 2 ) along the green-circle line. 
By direct inspection of the joint variance, this Aepr < 1 
already suggests that the two modes are entangled. 

Instead of extracting E w from a particular point in 
Fig. [2][e), we can reduce uncertainty by calculating 


by using the quadrature measurements made over all 
(01,02)- By fitting all the measured variances to the 
single squeezer model represented in Fig. [TJ we calcu¬ 
late the covariance matrix from extracted model param¬ 
eters (see Appendix ©. Consistent with the measured 
variances (Fig. [2j), the elements of the covariance ma¬ 
trix show both modest squeezing within each mode, and 
stronger intermode correlations (Fig. [3]). From the co- 
variance matrix, we directly calculate the entanglement 
witness E’w = —0.263 ±0.001 with a = 1.11. To quantify 
the degree of entanglement, we also calculate the nega¬ 
tivity N = 0.0824 ±0.0004. The uncertainties of Aw and 
N are estimated from a parametric bootstrap method. 
In the bootstrap procedure, 20 simulated data sets are 
generated from the measured covariance matrix. We es¬ 
timated statistical uncertainties as equal to the standard 
deviation of the 20 estimates Aw and N found by ana¬ 
lyzing the simulated data sets in the same manner used 
for analyzing real data. Furthermore, to demonstrate 
the repeatability and stability of the entanglement gen¬ 
eration, we perform 100 trials of the experiment. (The 
quoted Aw and N are one typical result from 100 tri¬ 
als of the experiment.) The distribution of Aw and N 
calculated from 100 trials appear to be Gaussian with a 
mean ± standard deviation of Aw = —0.264 ± 0.002 and 
N = 0.0820 ± 0.0007. 

Because we acquire a large data set in a short time, the 
statistical uncertainties of Aw and N are low compared 
to the systematic errors. We investigate systematic errors 
by analyzing experimental data with a general Gaussian 
model (see Appendix E]). Whereas the single squeezer 
model restricts possible estimated covariance matrices to 
those generated by a single mode squeezer, beam splitter, 
and loss, the Gaussian model allows any two mode state 
described by a Gaussian Wigner function. Using this 
Gaussian model, we found that Aw = —0.297±0.002 and 
N = 0.0921 ± 0.0004, where the uncertainties are based 
on a similar parametric bootstrap method to that de¬ 
scribed above. The discrepancy between the parameters 
estimated by the single squeezer model and the Gaus¬ 
sian model reveals systematic error that is significantly 
larger than any statistical uncertainty and uncertainty in 
the variances of the calibration state. Considering this 
systematic error, we conclude that the state in our ex¬ 
periment has Fw = —0.263^0.036 an d N = 0.0824^o.ooo4> 
where the uncertainties include the estimates from both 
models. 

The observed negativity is small compared to A = 0.61 
and N = 0.55 from [6] and Q, respectively, but those 
negativities are inferences of the negativity in the ab¬ 
sence of measurement inefficiencies. However, reference 
@ effectively combines two squeezed states on a hybrid, 
therefore, it is possible to generate states with larger neg¬ 
ativity using that method than the method demonstrated 
here and in [7]. We state the negativity without cor¬ 
recting for measurement inefficiencies. In comparison to 
the noisy amplification used in [tJ and the approximately 
14% efficiency achieved in [6], our apparatus achieves 
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measurement efficiencies at least 771 = 26 db 0 . 1 % and 
r ]2 = 41 db 0.2% (To quote these efficiencies, we assume 
t = 0.51 based on the calibrated measurements of the hy¬ 
brid [25]). Both this work and j!J benefit from quantum- 
efficient preamplifiers. The efficiencies presented here are 
higher than in [o] simply because we operate our pream¬ 
plifiers (VERs) with higher gains further reducing the 
apparent noise added by the following HEMT amplifiers. 
The efficiencies include any noise in the squeezed state’s 
generation and any loss along the entire path from state 
generation to measurement, in other words, the quoted 
values are lower bounds of the detection efficiencies. Be¬ 
cause we know only lower bounds on the detection ef¬ 
ficiencies we are unable to estimate the negativity (or 
other properties) of the generated state in the absence of 
measurement inefficiency. 

In conclusion, we demonstrate a two-channel, single¬ 
quadrature quantum measurement apparatus in the mi¬ 
crowave regime, where each channel of the apparatus uses 
a JPA as its first stage amplifier. A two-mode entan¬ 
gled state, which is generated by combining a squeezed 
state and vacuum on a microwave hybrid, is measured 
with improved efficiency and independent choices of each 
mode’s measured quadratures. Entanglement is demon¬ 
strated by showing that the two-mode state violates the 
separability criterion. Our integration of JPAs for both 
the preparation and measurement of an entangled state 
is a substantial addition to the toolbox for manipulating 
continuous variable quantum states of microwave modes. 
The measurement scheme is promising for demonstrating 
protocols exploiting entanglement. 


Appendix A: Vacuum calibration 

When we demonstrate the two-mode entanglement ex¬ 
periment, we feed the entangler with two input states 
emitted from the two 50 Cl terminations thermally an¬ 
chored to the cryostat. In additional to these two in¬ 
put states, more thermal (nearly vacuum) modes dilute 
the entangled state through the losses of the commercial 
microwave components, such as directional couplers and 
circulators. The temperatures of the input states and the 
loss modes are assumed to be equal to T in . To measure 

I 


Ti n , we inject a series of known thermal states into the 
two measurement chains without pumping SQ and calcu¬ 
late the variances Var(Vi) and Var(V 2 ) of the two output 
measurements for each thermal input. We then fit the 
model 


Var(V) = G 



~\~ A (Tp) , 


(Al) 


to the measured variances, where G and A(Tp) are the 
power gain and the added noise of the measurement 
chain. In this model, Ti n = ^/T| db T e 2 represents the in¬ 
put states temperature. Tp is the cryostat temperature, 
and we add one parameter T e to allow for the possibility 
that the terminations equilibrate at a higher tempera¬ 
ture than Tp. We also include a temperature dependent 
added noise A(Tp) = A 0 + A 2 TI in the model. The origin 
of the temperature dependent added noise in our JPAs is 
still under investigation, but it seems to be caused by the 
presence of resistive filters in the on-chip bias lines. From 
the fit (Fig. |4j), we extract 0 < T e < 16.1 mK. Because 
the entanglement generation is performed with the cryo¬ 
stat below 25 mK, we have T- xn <29.7 mK, which means 
the input variance is indistinguishable from vacuum in 
our experiments. 

The experiment’s data set contains amplified quadra¬ 
ture measurements Id,off (measured by bypassing SQ) of 
the vacuum state and quadrature measurements n of 
the entangled state for modes i = 1,2. Those variances, 
measured by digitizing the amplified voltages at room 
temperature, must be calibrated in units of quantum vac¬ 
uum at the input of the experiment. To calibrate the 
variances of the measurements in units of vacuum, we 
first normalize the quadrature measurements of the two 
entangled modes by the variances of the vacuum states 
to get U\(Qi) and 1 / 2 ( 02 ) 


UM 


Vl,on( 0 l) 


u 2 (9 2 ) 


V,on(fl 2 ) 

VMM' 


Furthermore, we derive the equations for variances of 
of U 2 (9 2 ), and of Ui(9i)± U 2 {0 2 ) from the single¬ 
squeezer model 


Var \Ui(0\)\ 

Var [U 2 (9 2 )\ 
Var \Ui(9i) ± U 2 (9 2 )] 
±^/gi 92 


9i 


92 


1 (s — 1 1 — 1 

1 + -- a + -- a cos(2#i + 20i) 


2 s 2 s 


2 1 s 2 - 1 

/? + -— 7 —/ 3 cos( 2 # 2 + 2 ^ 2 ) 


Var(17i) + Var (U 2 ) 

—-— \Ja/3 cos(#i + 0 2 + (j)i + <p 2 ) + -—-—— ya/3 cos (#2 — 9\ -\- (j) 2 — <f>i) ]. 

I- 


(A 2 ) 

(A3) 

(A4) 


The parameters a and j3 combine the hybrid power cou- pling coefficient t and the two measurement efficiencies 
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T CJyo (K) 
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FIG. 4: Thermal sweep experiment. The measured single-quadrature variances of the thermal sweep experiment (blue circles) 
and ht to Eq. (I All) (red solid line) for channel 1 (a) and channel 2 (b) are plotted. The residuals of the two fits are plotted 
in (c) and (d) where the error bars show the standard deviation of five independent measurements made at each temperature 
point. 


r]i and t] 2 

a = tr]i and /?=(! — t)r} 2 - 


Appendix B: Covariance matrix, the entanglement 
witness, and the negativity 


We do not assume that the VERs are perfectly linear, 
but rather introduce two parameters gi and g 2 to model 
the changes in VERs’ gains when bypassing SQ: 

9l,on 92, on 

9i = -, 9 2 =-, 

9 l,off 9 2,off 

where gi, on and gi, 0 ff are quadrature power gains of 
mode 1 for SQ is operated or bypassed and the same 
for mode 2. Finally, from the joint fit of the model equa¬ 
tions (lX2l)-(T£4l) to the measured Var(E/i), Var([/ 2 ), and 
Var(F r i =b U 2 ) [Fig. [5](b)— (d)], we extract s = 5.41 zb 0.03, 
a = 0.1304 zb 0.0007, f3 = 0.202 zb 0.001, 0i = -1.070 zb 
.0.002, 02 = -0.176 zb 0.001, g x = -1.70 zb 0.07 %, and 
g 2 = 2.04 zb 0.08 %. We then are able to use < 71 , g 2 , and 
Ti n to calibrate the quadrature measurements in units of 
vacuum: 


WxiOx) 


U\(0\) CJ 

V9i 0.5’ 


w 2 (e 2 ) = 


U 2 (02) a 

v. 0.5’ 


where a = (1/2) coth(ft/ s /2fcBT in ). 


To compute the covariance matrix predicted by the 
single-squeezer model, we examine the effects of each 
of the linear optical transformations (including squeez¬ 
ing) shown in Fig. 1 of the main text. Each of 
the transformations evolves the quadrature vector Z = 
{Xl, Yi,X 2 ,T 2 } t to MZ. The matrix M that describes 
the transformation is a real 4x4 matrix in the symplec- 
tic group Sp(4,M) [HI. The quadratures of the trans¬ 
formed state will have the covariance matrix MSM T , 
where E was the covariance matrix of the original state. 
The squeezing of mode 1 by amount s is described by 


S(s) 


A 

0 

0 

0 


0 0 0 

v^oo 
0 1 0 
0 0 1 
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10,000 samples, 1 ms 10,000 samples, 1 ms 10,000 samples, 1 ms 


(b) 



Time (ms) 

FIG. 5: Joint fit of the model equations to the measured variances, (a) The cartoon diagram represents the timing of the 
quadrature measurements. Each measurement is a one second long time trace with each channel sampled at 10 MSamples/s. 
Because the pumps of VER1 and VER2 are detuned from the SQ pump by 1 kHz and 50 kHz, respectively, the two VERs 
amplify the measured states at 10,000 different quadrature phase combinations every 1 ms (one record). In a 1 s long time trace, 
we thus have 1000 records, or realizations at each quadrature phase pair. The variance at each quadrature phase pair is then 
calculated from the corresponding data points in each record. The measured variances of (b) of (c) ^ 2 (^ 2 ), and of (d) 

1/2 [Ui(0i) + U2(02)] are plotted versus time in each measurement record (blue dot) along with the joint fit of Eq. (IA2D — (IA4D 
to all of the measured variances (red solid line). The model and data show good agreement with each other and we observe 
moderate squeezing below vacuum fluctuation (0.5) in the variance of each mode [(b) and (c)], but more squeezing in the joint 
measurement (d). 


A beam splitter of transmissivity t is described by 


Phase shifting mode 1 by p 1 and mode 2 by p 2 is de¬ 
scribed by 



r vi 

0 

— \/l — t 

0 


" cosOi) 

sin(0i) 

0 

0 ' 

B(t) = 

0 

\/l — t 

Vt 

0 

0 

Vt 

— y/l — t 
0 

P(4>i,h) = 

— sin(</>i) 
0 

cos(0i) 

0 

0 

cos((p 2 ) 

0 

sin(0 2 ) 


0 

y/l-t 

0 

Vi 


0 

0 

-sin (> 2 ) 

cos (<p 2 ). 
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Beginning with an initial vacuum state, Ho = 1/2, 
(with I being the identity matrix), we compute 
the state just prior to measurement to be E p = 

To account 

for photon loss we append two ancilla modes, which are 
coupled to modes 1 and 2 with beam splitters B(rji) and 
B{r] 2 ) respectively. After this coupling, the ancilla modes 
are discarded. This procedure transforms E p into 

H ss = HTi p + H- } 

I 


where H = diag (771 , 771 , 772 , 772 ) and H = diag(l — 771 ,1 — 
V 1,1 —772,1 —%)• 


Given a covariance matrix E, we compute the value of 
the entanglement witness Pw (shown in Eq. (j2j) of the 
main text) by 


E w = min 

Oi,02,cl 


a x T P(e 1 , 0 2 )XP(e 1 ,8 2 ) T a x + a y T P(6 1 ,8 2 )j:P{e 1 , e 2 ) T a y - a 2 + 


where a x T = (|a|, 0,1/a, 0) and a y T = (0, |a|, 0, —1/a). 
In fact, P w ’s only dependence on 0 1 and 62 appears in 
the form 6 \ — 62 , so it is only necessary to minimize over 
one phase. 

To compute the entanglement negativity from a co- 
variance matrix, we follow the treatment given in poj . 
Through the application of linear optical devices and 
squeezing (the symplectic transformations), any 2 mode 
covariance matrix E can be transformed into the co- 
variance matrix of a thermal state, which has the form 
diag(zq, zq, ^ 2 , ^ 2 )- v\ and V 2 are called the symplectic 
eigenvalues of E. According to the Heisenberg Uncer¬ 
tainty Principle, v\ and ^2 > 1/2. If the quantum state 
p with covariance matrix E is separable, then covariance 
matrix E of the partial transpose of p will also have sym¬ 
plectic eigenvalues v\ and > 1/2. If v\ or z/ 2 < 1/2, p 
must be entangled. To compute the symplectic eigenval¬ 
ues of E, we divide E into 2x2 blocks: 


E = 


a r 
r T b. 


The quantities |E| and A(E) = |A| + \B\ + 2|T| are invari¬ 
ant under the symplectic transformations. (| • | denotes 
the determinant.) From them we calculate the symplec¬ 
tic eigenvalues: 


Vi = \Jl (a(E) ± i/[A(£)] 2 — 4|5^j), 

where we use i = 1 for the — case and i — 2 for the 
+ case. The partial transposition of p has the effect of 
reversing the sign of the Y quadrature of the transposed 
mode, so that A(E) = \A\ + \B\ — 2|T|, but |E| = |E|. 
The symplectic eigenvalues of E are 


Vi 



(a(E)±0A(E)]2 



If v\ < 1/2, p is an entangled state. Finally, the negativ¬ 
ity is 


N = max (ft. ——. 

\ 2z/ i / 

Note that our formula for N is slightly different from that 
in [ 20 I , because that paper uses the convention that the 
variance of the vacuum state is 1, whereas we use vacuum 
variance of 1/2. 


Appendix C: Gaussian state estimation 

As a check for systematic errors we implemented a sec¬ 
ond method to estimate the quantum state produced in 
this experiment. We call this method “Gaussian state 
estimation”, because its range is all two-mode states 
that have Gaussian Wigner functions. As inputs it ac¬ 
cepts the calibrated quadrature measurements Wi( 0 i) 
for i = 1,2 and returns the Gaussian state’s vector of 
quadrature expected values p = (fei, p y i, p X 2 , l^y 2 ) T = 
E[(Xi, Yi, X 2 , Y 2 ) T ] (E[x] is the expectation value of 
x.) and the covariance matrix E whose elements 
are the covariances between of the random variables 
(Xi, Yi, X 2 , U 2 ) t . We will label the elements of E with 
subscripts indicating the quadrature variable and mode: 
E a i,bj, where a and b G {x, y} and i and j G {1,2}. Note 
that the Gaussian state estimation requires the use of 
the squeezer model described in Appendix lAl to produce 
correctly calibrated quadrature measurements. The data 
set consists of the n instances of the quadruplet con¬ 
taining two phases and two quadrature measurements: 
{(^ 1^5 w[ k \ 9^\ W 2 ^)\k = l,...,n} in which k labels 
the measurement instance. Although in the experiment 
Ox and 62 are scanned continuously and quadratures are 
measured at regular intervals, in this section we treat the 
phases as random variables uniformly distributed over 
[ 0 , 27r) with probability distribution P(0i) = l/(27r). 
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Consider the expected value 


E[Wi cos 6 i\ = E[(Xi cos Oi + Y{ sin Oi) cos Oi ]. 


Solving these three equations for the elements of E and 
applying the Law of Large Numbers gives us the esti¬ 
mates 


Because Oi is independent of Xi and 


E [Wi cos 0i\ = E [/jb X i cos 1 2 3 4 5 6 7 Oi + iiyi sin Oi cos Oi] 

27T 

{^/ixi cos 2 Oi + fiyi sin Oi cos Oi) P(0)d0 



E xi,xi = 3 Wf cos 2 Oi - Wf sin 2 Oi - ii x 2 
E yii yi = 3W 2 sin 2 Oi — Wf cos 2 Oi — p^ 2 
Yxi^yi = 4W^ COS Oi sin Oi l^xil^yi m 


Estimates for the cross-mode elements of E are given by 


By the Law of Large Numbers, we can estimate fi x i with 

1 n 

/Si = 2 Wi cos Oi = -Z W[ k) cos 6 {k) , 

11 k=1 

where we use the hat to denote the estimate of a param¬ 
eter and the overline to denote the sample mean. Simi¬ 
larly, we can estimate fi y i with jj^i = 2Wi sin Oi. Apply¬ 
ing this treatment to both modes gives us 
To estimate E, consider the expected value 

E [Wf cos 2 Oi] = E [(Xi cos Oi + Y t sin#;) 2 cos 2 0*]. 


Using the independence of 0; from X t and Pi and the 
uniformity of P(6i ), we obtain 


E[W 2 cos 2 9i] = E 


1(3 X? + Y?) 


g(3 






^xl,x2 = 4 VEi COS O1W2 COS O2 ~ T^xUPci 

E^i iV2 = 4Wi cos O 1 W 2 sin O 2 
E y i, s2 = 4LEi sin6>iVE 2 cos6> 2 - 
E y i,y 2 = dWi smOiW 2 sin0 2 - MyT/^ 2 - 

Thus, we estimate the expected values fi and covari¬ 
ance matrix E by computing sample means of simple 
functions of the quadrature and phase measurements. 
Because the computation is so simple, our method is 
well suited for large data sets. The maximum likelihood 
method developed by Rehacek and co-authors might give 
lower statistical uncertainty at the cost of greater com¬ 
putation time [26]. Our method does not impose a con¬ 
straint on E to enforce the Heisenberg Uncertainty Prin¬ 
ciple m. Although it is possible for the estimate to 
be unphysical, all states found in our analysis are valid 
quantum states. 


One can similarly show that 

E[fU 2 sin 2 Oi] = 2 xi + Sit, xi + 3fi 2 yi + 3T, yi , yi ), 
and 

E[W; cos Oi sin Oi] = ~X^l^ x il^yi T E x i,yi). 
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